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^ ; Abstract 

' This paper deals with dependence across marginahy exponentially distributed 

arrival times, such as default times in financial modeling or inter- failure times 

' in reliability theory. We explore the relationship between dependence and the 
Ph ! 

• , possibility to sample final multivariate survival in a long time-interval as a se- 

■ quence of iterations of local multivariate survivals along a partition of the total 
time interval. We find that this is possible under a form of multivariate lack of 

. memory that is linked to a property of the survival times copula. This prop- 

■ erty defines a "self-chaining-copula" , and we show that this coincides with the 
I extreme value copulas characterization. The self-chaining condition is satisfied 
. by the Gumbel-Hougaard copula, a full characterization of self chaining copu- 

o 



las in the Archimedean family, and by the Marshall-Olkin copula. The result 
I has important practical implications for consistent single-step and multi-step 

simulation of multivariate arrival times in a way that does not destroy depen- 
^ I dency through iterations, as happens when inconsistently iterating a Gaussian 

c3 , copula. 
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1 Introduction and Motivation 

This paper deals with dependence across marginally-exponential arrival times, for ex- 
ample default times in financial modeling or inter-failure times in reliability theory. 
Dependence can be modeled both through a multivariate distribution function consis- 
tent with the exponential marginal distributions, or by joining the given exponential 
marginal distributions through a copula function. After introducing a multivariate 
lack of memory property, we first present two multivariate exponential distributions 
that are known in the literature and satisfy such property, and then move to copula 
functions that are consistent with these distributions and with the lack of memory 
property in particular. The lack of memory property is important because it allows 
for consistency between a unique survival sampling at a single large time interval and 
iterated survival samplings at smaller sub-intervals forming a partition of the whole 
interval. Since in the industry it is important to have the possibility to decompose 
a final survival simulation into survival subsimulations, this is also a relevant prac- 
tical problem. The relevance stems from the need to make the survival simulation 
time step consistent with the common time step for all other underlying variables 
or risk factors (e.g. equity prices, interest rates, temperatures, etc). We provide 
a characterization of "self- chaining" copula functions satisfying the lack of memory 
property. We consider as fundamental examples the Marshall-Olkin copula and the 
Gumbel-Hougaard Copula. In the family of Archimedean copulas, we provide a full 
characterization, showing that the only lack-of-memory and consistent single-step 
and multi-step Archimedean copula is the Gumbel-Hougaard copula. 

The paper is organized as follows. Section [2] introduces the common period lack 
of memory (CPLM) property that makes consistency between single-step and multi- 
step sampling possible. The following Section [3] introduces multivariate exponential 
distributions that satisfy such lack of memory property. This is done directly at multi- 
variate distribution level without splitting marginal and dependence information, i.e. 
without resorting to copula functions. We resort to copula functions in the following, 
so that Section H] presents a quick introduction to copula functions, whereas Section 
[5] introduces the self-chaining property that translates the lack of memory and con- 
sistency above into copula language. Section |6] goes into more detail and derives the 
characterization of the common period lack of memory property in terms of copula 
properties, proving that this is achieved by the self-chaining property. A first char- 
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acterisation of self-chaining copulas in terms of homogeneity or PDEs is presented. 
Section [7] fully characterises self chaining copulas in the archimedean family as the 
Gumbel-Hougaard copula. Section [8] explains that self-chaining and extreme-value 
copulas are the same, and by summarizing the Pickands function characterization 
of Extreme Value Copulas provides a precise characterization also of self-chaining 
copulas. Finally, Section |9] concludes the paper. 

2 A multivariate form of lack of memory: CPLM 

Consider event times (typically inter-arrival times, inter-failure times or for example 
defaults) for n entities, and let the event time for entity i be denoted by r^. These 
event times can for example be default times of a number of possibly related firms or 
first-failure times for a number of entities. 

Denote by Fi the cumulative distribution function for r^, and by Gj = 1 — the 
related survival function. Suppose now that we look at a single period vs multi-period 
context. Consider a final period NT obtained by adding a partition of single period 
intervals of length T, hence periods [0, T), [T, 2T),...,[{N - 1)T, NT]. 

Assume that the single event times do not have memory, in that 



for any i and any < S < T (and hence are exponentially distributed and Gi{NT) = 
Gi{T)^). An important practical question is the following. Is it possible to iterate a 
simulation of the joint survival of arrival times always in the same way in all sub in- 
tervals [0, T), [T, 2T), . . . , [(A^ — 1)T, NT], and also in the same way as we simulate 
joint survival of arrival times in a single sampling run in [0,A^T]? The answer is 
affirmative under a multivariate form of lack of memory, namely (for the case n = 2 
for simplicity) 



P(ri > jT, r2 > jT\n > hT, > hT) = P(ri > (j - h)T, > (j - h)T) (2) 



with j > h integers. In this paper we will try and characterise this lack of memory 
property, which we term "common periods lack of memory" (CPLM), in accordance 
with the following remark. 

Remark 2.1. (Lack of memory at common levels for both random times) 

It should be noted, importantly, that we do not adopt the most general definition of 
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bivariate lack of memory, namely 

P(ri > jT, T2 > kT\n > hT, > iT) = P(ri > (j - h)T, T2 > (k - i)T) (3) 

for ma.x{i,h) < min{j,k). This is because we plan to apply our result to joint sim- 
ulations of the multivariate vector of arrival times in common intervals. Besides, 
Condition in Eq. ^ would be too strong, implying a trivial case of lack of mem- 
ory, namely independence of the two exponential random variables ti and T2, see for 
example 12^ . 

3 Multivariate exponential distributions and CPLM 

We now present two possible multivariate exponential distributions with exponential 
marginals and satisfying the Common Period Lack of Memory (CPLM) property. 

Consider the bivariate case for simplicity, involving a bivariate vector of univari- 
ate exponential arrival times ti,T2. Leaving aside the general metod of introducing 
a copula to connect marginally exponential random variables for the time being, it 
is well known that there are several bivariate exponential distributions that could be 
used to model the bivariate random vectors consistently with the univariate exponen- 
tial distributions. One of the most utilized bivariate exponential distributions is the 
Marshall-Olkin bivariate distribution (see [29]). This distribution generalizes the lack 
of memory property of univariate exponential distributions, satisfying in particular 
our CPLM definition in Eq ([2]) above. The Marshall Olkin bivariate distribution, 
however, features a singular component and admits a strictly positive probability 
that Ti = T2. When one excludes perfectly simultaneous ti and T2, one may resort for 
example to a bivariate exponential distribution among the three proposed by Gumbel 
in [H]. The first bivariate distributions in [ll] satisfies an alternative definition of 
lack of memory, also known as bivariate remaining life constancy, see [21]. However, 
this first distribution proposed by Gumbel can only describe negative dependence and 
in a limited range. The second bivariate exponential in pJ:J only describes a range 
of dependence [—1/4,1/4] for correlation, and as such is of limited scope. This is 
why we resort to the third bivariate exponential distribution only briefly introduced 
in [H]. See also [25] and [21]. 

The joint survival function for this bivariate exponential is, for simple univariate 
exponential marginals with positive intensities A's, and for dependence parameter 
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Q(ri > xuT2 > X2) := G{xi,X2) = exp(-((AiXi)' + (Aa^a)')'/'). (4) 

Notice that, indeed, the marginal distributions are exponential random variables 
with mean respectively 1/Ai and 1/X2- We will set Ai to a constant deterministic 
intensity, and A2 to a constant intensity as well. It is straightforward to check that 
this third Gumbel multivariate exponential distribution satisfies the lack of memory 
property in Eq ([2]) (i.e. CPLM, but not the more general lack of memory property 
in EqEl). 

We notice that Kendall's tau for this distribution, which is a good measure of 
dependence (invariant for invertible increasing transformations), is 

T^{G) = 1 - 1/e. (5) 

This confirms that 6 = 1 characterizes the independence case, whereas — ?■ 00 char- 
acterizes the co-monotonic case. This bivariate exponential distribution can therefore 
describe the whole range of positive dependence. 

We can also notice that A's are pure marginal parameters, whereas 6' is a pure 
dependence parameter. This is a bivariate distribution allowing for tail dependence, 
and, differently from Marshall Olkin's, does not have a singular component, so that 
there is zero probability that ri = T2. 

It should be noted that in a number of applications, and especially in financial 
modeling of default times, the assumption that A's are deterministic is not realistic 
and its negative features have been highlighted for example in [5]. 

The Marshall- Olkin and Gumbel multivariate exponential distributions satisfy the 
CPLM property and thus lend themselves to a consistent single-step and multi-step 
simulation and sampling of survival. However, to better characterize dependency 
that is consistent with multivariate lack of memory we standardize away marginal 
distributions and focus on copula functions as follows. 

4 Copula Functions 

Consider a random vector X = (Xi, X„), and suppose that we wish to analyze the 
dependence between its components. The law of X is then characterized by its joint 
cumulative distribution function (CDF) {xi, . . . ,Xn) ^ IP(Xi < xi,...,Xn < Xn)- 
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However, this function mixes information on the dependence between the different 
components of the vector with information on the distribution of the single com- 
ponents themselves. Copula functions have been introduced in order to allow for 
a separation between the marginal CDF and the dependence structure. The for- 
mer concerns single components, taken one at the time, and is given by the CDF's 
Fi{x) := < x), i = l,...,n, which we assume to be continuous. The latter 

is entirely represented by the copula function we introduce now. It is well known 
that Ui = Fi(Xi), f/„ = Fn{Xn) are uniformly distributed random variables on 
[0,1]. The joint cumulative distribution function of {Ui, ...,Un), that we denote by 
C{ui, ...,Un) = P(f/i < Ui, ...,Un < Un), is Called the copula function of (Xi, ...,X„) 
and satisfies: 

P(Xi < Xi, X„ < Xn) = C(P(Xi < Xi), P(X„ < Xn)). (6) 

A copula function has the following three properties: 

C{Ui, ..,Mi_i,0,Mi+i, ..,Un) = (7) 

C(l,..,l,Ufc,l,..,l) = Mfc (8) 
dui...unC is a positive measure in the sense of Schwartz distributions. (9) 

Condition ([9]) means concretely that for any hypercube 

H = [ai, 61] X ... X [a„, 6„] C [0, 1]" we have F[{Ui, .., f/„) eH]>0. 

Conversely, one can show that any function that satisfies the above three condi- 
tions (iTl [HI El) can be viewed as the joint CDF of a vector of uniform variables on 
[0, 1] and is thus a copula, see for example [16], [31], [32] and [M]. In this spirit, in 
the following, the expression "simulating a copula C" will denote the simulation of a 
random vector of uniform variables {Ui, .., Un) on [0, 1] whose joint CDF is C . This 
can be occasionally referred to as "sampling the copula". 

Among the different ways to define specific copula functions, there are the follow- 
ing two. The first one consists in seeking functions C satisfying the three above prop- 
erties. Archimedean copulas are an example of this approach. Indeed, Archimedean 
copulas stem from the remark that if is a convex strictly decreasing function such 
that <y9(l) = and with suitable additional properties (see Theorem 17. II below) . then 
C(mi, ..,Mn) = </''~^'(v^(wi) + .. + 'f{un)) has the above three properties and is thus a 
copula. The function (p is called the generator of the copula. The Archimedean copula 
is said to be strict if the generator satisfies lim^i_j,o+ ^p{u) = 00, and for these copulas 
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the pseudo-inverse yjt"^' can be replaced by the inverse y?"^ and the generator is said 
itself to be strict. Examples of Archimedean copulas are the Gumbel-Hougaard, Joe, 
Clayton, Frank copulas, see Nelsen (1999) and Joe (1997). 

The second method consists in working directly with joint CDF's F{xi, ...,Xn) and 
the related marginal CDF's Fj, providing the copula F{F^^{ui), F~^(m„)). Even 
if this method does not always lead to analytically tractable copulas, it can provide 
us with copulas that are easy to simulate. Indeed, the main example of this kind of 
construction is the fundamental family of Gaussian copulas. Student-t copulas are 
also defined this way. 

A Gaussian copula is defined as the copula of a joint Gaussian random vector 
X with standard Gaussian marginals and correlation matrix p, and is thus given by 
Np{N~^{ui), N~^{un)) where N is the CDF of a scalar standard Gaussian variable 
and Np is the joint CDF of X. 

In recent years, copula functions have received a great deal of attention, see for 
example the papers of Genz and Bretz (2002), Hiirlimann (2003), Juri and Wiithrich 
(2002), Wei and Hu (2002), Alfonsi and Brigo (2005), and the books of Joe (1997) 
and Nelsen (1999). For financial and insurance applications, recent applications on 
copulas include for example Li (2000), Bouye et al. (2000), Cherubini et al. (2002), 
Jouanin et al. (2001), Klugman and Parsa (1999). For an excellent extensive overview 
see Embrechts et al. (2003), and for applications to finance see [8]. 

5 Multivariate Sampling: Single vs Multi step cop- 
ulas 

Consider the joint survival or no-event (no-failure) probability function at ti, . . . ,tn, 

F(ri >h,...,Tn>tn)= P(G'i(ri) < G'i(ti), . . . , Gn{Tn) < 

= < Gi(ti), . . . , f/„ < GM) =: C{uu ...,Un) 

where we termed C the survival times copula and we set Ui = Gi{Ti), uniform random 
variable, and Ui = Gi(ti). 

Suppose now that we look again at a single period vs multi-period context. We 
still enforce univariate lack of memory as in Eq ([T]). In this paper we will try and 
characterise the multivariate lack of memory property CPLM of Eq ([2]) in terms of 
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copula properties. We will see that if CPLM holds, we can sample the event where all 
arrival times exceed time NT either at the final time NT directly, which is consistent 
with a survival times copula 

CK, . . . u, = G'i(7VT), ...,Un^ Gn{NT), 

or by iterating steps of survival in each [{i — 1)T, iT), i — 1, . . . ,n. We will see 
that under CPLM the second procedure produces the copula 

Clearly, to have consistency the two copulas have to coincide, namely 

{C{uY'',...yJ'')f = C{u,,...,Un) 

When this happens for a given copula function, we call the copula a "self-chaining 
copula" , and this is consistent with common period multivariate lack of memory. 

Now the problem is that the industry wrongly applies the above procedure or 
related procedures even when the copula is not self-chaining and when CPLM does 
not hold. Typically, the above procedure is applied with a Gaussian copula C. This 
is not consistent and leads to a difference between the one-shot sampling at NT and 
the iterated sampling, to the point that the latter kills dependence for large N. 

This is such an important point that it is worth clarifying it with an example. 

Example 5.1. (Iterating a Gaussian copula kills dependence) Consider two 
exponentially distributed default times connected by a Gaussian copula with depen- 
dence parameter p and with constant intensities Ai and X2, namely 

n = - ln(l - $(Xi))/Ai, T2 = - ln(l - $(X2))/A2, 

where [Xi , X2] is a bivariate Gaussian random vector with standard gaussian marginals 
and with correlation parameter p. 

Assume Ai = A2 = 0.02, p = 0.9, and consider the following two procedures: 

• Sample directly the event Ti > lOOy Cl T2 > lOOy. The probability of this event, 
based on a simulation with one million scenarios, is 



0.097 ±0.0003 (0.0969) 
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where we included the simulation standard error. The value within brackets is 
the value computed by resorting to a double numerical integration routine to 
compute the Gaussian copula. 

• Iterate the check AjTi > ly (1 Ajr2 > ly 100 times, where [A^ri, Air2], i = 
1,2, ...,100 are independent copies of [ti,T2] to be used to check default in 
every year. We count the scenarios along which there is always joint survival 
up to 100 years and divide by the total number of scenarios. 

This yields a survival probability of 

0.057 ±0.0003 (0.0557) 

The value within brackets is the value computed by resorting to a double numeri- 
cal integration routine to compute the Gaussian copula for ly, thus determining 
the ly survival probability, and then using the 100th power of this number for 
the survival over 100 years. 

As we can see, there is a quite relevant difference, in this case, between simulating 
joint survival one-shot for 100 years, and iterating 100 times a ly survival event as 
if the multivariate arrival times satisfied lack of memory. In the latter case the joint 
survival probability (and hence dependence) is much smaller. This is a confirmation 
that iteration, when lack of memory is not satisfied, kills dependency. This happens, 
in particular, with the Gaussian copula. 

In this paper we will characterize copula functions leading to CPLM and hence 
to the possibility to have consistent sampling and subsampling, or single-step and 
multi-step simulations which do not destroy dependence through iteration. This is of 
great practical relevance. 

6 Copula characterization of multivariate Lack of 
Memory 

The main result of this section is the following 

Proposition 6.1. Assume single event times ti and T2 satisfy each the lack of mem- 
ory property (0]). Lack of Memory extends to the bivariate default time via the 
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CPLM condition in Eq ^ if and only if the survival or no-event time copula satisfies 
C{ui,U2) = C{ui,U2)^ (self-chaining copula) (10) 
for all ui,U2 in [0, 1] and integer positive k. 

Proof. Assume lack of memory holds in the bivariate case too. Compute 

P(ri > 2T, T2 > 2T) = P(ri > 2T, > 2T|ri > T, ra > T)P(ri > T, Ts > T) = 
= P(ri > 2T - T, Ta > 2T - T)P(ri > T, rs > T) = P(ri > T, ra > T)^ 

Now lack of memory for single default times implies easily Gi{2T) = GilT)"^, so that 
by taking G"s on both sides of the events whose probabilities are computed, from the 
previous equation we have 

P(f/i < G,{T)\ U2 < G2{Tf) = P(f/i < G'i(T), U2 < G2{T)f 

or actually ffTOl) with k = 2. Iterating shows ffTOl) holds for all k. 

Viceversa, if ([TI]]) holds for all k, it is immediate to prove that lack of memory 
holds. □ 

Remark 6.2. (Extension to real exponents) 

The change of variables Ui^2 '■= 2 with a substitution in shows that the anal- 
ogous of nOj) in Ui^2 will hold with exponent 1/k. Combining the two results for 
different k yields easily that fT^) will hold for any rational positive exponent replac- 
ing the positive integer k, and by extending the relationship by continuity to the whole 
set of positive real exponents we have [W\] holding for any positive real k. 

The proof also extends straightforwardly to the copula on n exponentially dis- 
tributed arrival times. 

We refer to Condition ( ITOj) . that in the n-dimensional case reads 

C(m^, • • • 5 ^n) = C{ui, U2,..., UnY, G (self-chainiug copula), (11) 

as to a "self-chaining" property of the copula. This property tells us that if we break 
the simulation of first arrival times up to a final time jT into simulation on intervals 
[1,T], [T, 2T],..., [(j — 1)T, jT] then both the overall copula in the one-shot case and 
the step by step copula in the multi-step case are the same copula and lead to the 
same final result. 

A simple characterization of the self-chaining condition guaranteeing common 
periods lack of memory for the multivariate exponential arrival-times is the following. 
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Proposition 6.3. Let L{vi,V2) '■= \ogC{e'"^,e^'^), i.e. the log-copula with exponen- 
tially expressed arguments, defined for vi,V2 such that C{e^'^,e'"^) > 0. We have 

Ciul 4) = C(mi, M2)^ ^ L{kv^, kv2) = kL{vuV2) (12) 

for positive real k and all Ui,U2 where C is strictly positive. 

This amounts to homogeneity of order 1 in the positive real k for the log copula 
with exponential arguments. The proof is immediate. 

It is easy to notice that, for known families of copulas, Property ( fT2l) is satisfied 
by the Marshall-Olkin copula and the Gumbel-Hougaard copula. 

Proposition 6.4. The Gumbel-Hougaard copula and the Marshall-Olkin copula are 
self- chaining. 

Problem 6.5. (^Characterization of Self Chaininig Copulas) Find a character- 
ization of self-chaining copulas, including Marshall-Olkin and Gumbel-Hougaard as a 
special case. In a way, this amounts to characterize homogeneous functions L (i.e. 
satisfying Eg [71|) whose exponential exp(L) satisfies Eqs. (0 0, 

The homogeneity condition may also be expressed as a PDF. Starting from 

L{tx,ty) = tL{x,y), 

and taking the total derivative wrt t on both sides one obtains 

L^{tx, ty)x + Ly{tx, ty)y = L(x, y) 

where Lj.{x,y) = dxL{x,y) and similarly for Ly. Set t = 1 to get 

L^{x,y)x + Ly{x,y)y = L{x,y). (13) 

Remembering that L(t>i,t>2) := \ogC{e^^ ,e'"^) one has 

Cu{u,v) u \og{u) + C^iu^v) V log(f) = C(-u,f) logC(M,f). 

Now the problem is: Characterizing all copula functions, i.e. all positive functions C 
with values in [0, 1] and satisfying Eqs (0 [3 that satisfy this PDE. 

Alternatively, one can work with the L PDE [T^) . with L taking values in (— oo, 0], 
and such that exp(L) satisfies Eqs (0 0, [3^. 
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In the strict Archimedean copulas context, with copulas associated with frailty 
distributions, we can solve the above problem and provide a full characterization of 
the self chaining property. 

Before solving the above problem in the archimedean context, it is worth noticing 
that the self-chaining Marshall-Olkin and Gumbel-Hougaard copulas are associated 
resprectively with the multivariate Marshall-Olkin and Gumbel (Eq. (^) multivariate 
exponential distributions mentioned in Section |3] when one standardizes away the 
marginals. Such distributions satisfy the CPLM property as we had already observed 
in the abovementioned Section. 

However, despite such consistency, the Marshall-Olkin and Gumbel-Hougaard 
copulas are far from being the most commonly used copula function. The most 
commonly used copula functions are the Gaussian copula and the Student-t copula. 

Gaussian Copula functions have been associated by part of the press to the fi- 
nancial crisis started in 2007. For a discussion on a number of misunderstandings 
associated with this opinion and on the real shortcomings of the way Gaussian copulas 
have been used for Collateralized Debt Obligations (CDOs) see [6]. 

Remark 6.6. (Iterating a Gaussian copula for multivariate arrival times 
sampling is inconsistent). It is worth saying immediately that the Gaussian copula 
is NOT self-chaining and therefore is not consistent with multivariate lack of memory. 
Iterating a Gaussian copula for joint arrival times simulation is thus wrong, and kills 
dependence in the long run, as we have seen in Example \5.1\ above. 

In the following we turn to characterization of self-chaining copulas in the Archimedean 
family. 

7 Characterizing Archimedean self-chaining copu- 
las 

To accomplish the above characterization within Archimedean copulas, we need to 
narrow a little the focus to Archimedean copulas whose generator is associated to a 
Frailty distribution. The key result in this respect is Kimberling's (1974) characteri- 
zation of Archimedean copulas with a generator working for every dimension: 
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Theorem 7.1. (Kimberling, 1974) Given a strict generator ip, 

the function {ui, . . . , Un) i->- (fi~^{(p{ui) + .. + (fi{un)) from [0, 1]" to [0, 1] is a copula 

for every n if and only if (fi~^ is completely monotone on [0, oo). 

In turn, the Bernstein- Widden Theorem (Widder 1946) states that complete 
monotonicity is equivalent to being the Laplace transform of a non-negative Ran- 
dom variable measure, that is called the frailty distribution. Therefore, given a strict 
generator (p, this generates a strict Archimedean copula for every dimension n if and 
only if the inverse generator 99"^ is the Laplace transform of a frailty distribution. A 
discussion on this characterization and some interesting cases where complete mono- 
tonicity does not hold are in McNeil and Neslehova (2007). 

Before we proceed further we also need some definitions and useful propositions. 
We follow Sato (1999) and for a probabihty measure // we define its Fourier transform 
with ji. 

Definition 7.2. A probability measure fj, is infinitely divisible if for each N G N 
there exists a probability measure /in such that 

A = (Aiv)^ 

Theorem 7.3. // an Archimedean copula is self chaining, then its frailty distribution 

is infinitely divisible. 

Proof For all z e [0, +00) set u = (f)'^{z/2) e (0, 1]. Now write 

0-i(^)Viv ^ 0-i(20(^,))Viv ^ 0-i(20(x,^)) = 0-i(20(0-^(z/2)V^)) 

Therefore, if we set (pN^iz) = 0-1(20(0-^(^/2)1/^)) we can write 

r\z) = rN\zr 

Since the inverse generator is the Laplace transform of the frailty, the same rela- 
tionship will hold for the characteristic function. Therefore the frailty is infinitely 
divisible. □ 

Definition 7.4. Let /i be an infinitely divisible probability measure or R. It is called 
strictly stable if for all a > 0, there is b > such that 



ji{zr^ji{bz). 
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It is called strictly semi-stable if for some a > with a 1, there is b > satisfying 
the above relationship. The stability index is the value a that solves 

a\b\-"^l 

Theorem 7.5. If an Archimedean copula is self chaining, then its frailty distribution 
is strictly semi-stable. 

Proof. Consider an Archimedean copula with generator (/?, and define a set if indices 
/. Then, for a vector u = {ui}i^i e [0, 1]^ the copula will be given by 

For the copula to have the self-chaining property we demand that for all £ > 

c {{uIUj) = c i{u,},^jY 

For each N e N we therefore compute the copula value at a point with N abscissas 
equal to u and the rest equal to one. Thus, the self-chaining property implies 

We now define an auxiliary function g and its inverse as 

g : [— oo, 0] [0, -|-oo] : x 0(expx) 
g-^ : [0,+oo] ^ [-00,0] : y log0-^(y) 

Taking logs of both sides of the self-chaining property allows us to write it in terms 
of g in the following way 

g-^ {Ng{ix)) = ig-' {Ng{x)) 

Thus, if we define a second auxiliary function Gn{x) = g^^{Ng{x)) the above rela- 
tionship becomes 

Gn{£x) = £Gn{x) V£ > 0, X < 

This implies that Gjv is a linear function without a constant. To see that, set x = — 1, 
z — —i and kN — Gn{—1). Then write 



Gn{z) = -zGn{-1) = -kNZ < 
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HavingestablishedtheformofGiv we substitute backwards, writing ^{Ng{z)) — 
—kj^z, and by the change of variable x — Ng{z) we write that 

g-\x) = -kng-\x/N) 

As g^^{x) = \ogip~^{x) we conclude that self-chaining implies a generator that satis- 
fies the relationship 

[^-\x)]-'/'^^^-\x/N) 

The inverse generator is the Laplace transform of the frailty distribution. The 
Fourier transform will be given of course by fi{z) = (j)~^{—iz). The relationship will 
then carry onwards to the characteristic function of the frailty. 

We have therefore shown that there exist epochs ajv = — I/Zcat, for which we set 
spans — ^/N, satisfying the relationship 

Fix)"'' = F{bNx) 

By definition we conclude that the frailty distribution is strictly semi-stable. We 
denote the stability index with a. □ 

Prom the definition of the stability index we can express the relationship between 
the spans and their corresponding epochs. In particular 

It is possible to characterize the frailty distribution even further, since the gen- 
erator of the frailty distribution is continuous. To do so we also need to define the 
measure operator for each r > 0, which transforms a measure // on IR as 

{Trn){B) = //(r-^S), for Borel B e B(R) 

In addition we provide a useful theorem (Sato, Thm 14.3) 

Theorem 7.6. Let be infinitely divisible with Levy measure v. Let < a < 2. 

1. II is a-semi-stable with b > 1 as a span iffu satisfies 



2. ji is a-stable iff the above relationship holds for all b > 0. 
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Theorem 7.7. If a strict Archimedean copula is self chaining then its frailty distri- 
bution is strictly stable. 

Proof. We have shown that the fraihy is strictly semi-stable, with spans = 
and associated epochs qn = —(1/-^)", for all iV G N. It is also straightforward 
to verify that the quantities b^ — N are also valid spans, with associated epochs 

Therefore, for any rational number b — M/N and any Borel B e B(]R) we can 
write 

inu) {B) = i^-j (Tm/at^) (B) = ^-j {TmT,/^v) (B) = u{B) 

Since the distribution is strictly semi-stable for all spans b = M/N, M, N & N, 
by continuity it is strictly semi-stable for all spans 6 > 0, and is therefore strictly 
stable. □ 

Finally, we can characterize continuous self-chaining copulas completely. 

Theorem 7.8. // a strict Archimedean copula is self chaining, then it is the Gumbel- 
Hougaard copula. 

Proof. Having established that the frailty is strictly stable with positive support we 
can write the characteristic function in closed form. This will in turn determine the 
generator. 

In general, stable distributions will have characteristic functions given by 

L{—iz) = exp |— c|^|° ^1 — i/9tan ^ sgn^;^ -|- i72;| , if a 7^ 1 

L{—\z) — exp |— c|z| ^1 + i/9— sgn 2;log -|- vyz^ , if a = 1 

for c > 0, /3 e [^1-1] and 7 G R. These three parameters, together with the stability 
index a determine the density. 

As the frailty has positive support the restrictions /5 = 1 and 7 = are im- 
posed. Under these constraints the Laplace transform, which will also be the frailty 
generator, will be defined only for a 7^ 1, and is given by 

L{z) — exp |— c|2;|" ^1 -|- tan | , for a 7^ 1 

Also, as the Laplace transform is decaying the further constraint 

ttq; 

1-Ftan — >0=^q;<1 
2 
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In this case the quantity 5 = c (l + tan ^) > is a constant, and the inverse 
generator takes the Gumbel-Hougaard form. □ 

8 Relationship with extreme value copulas 

We found out later, after proving the above result, that such result is already known 
in a different context. In particular, it is a result associated with max stable copulas 
(MSG) and Extreme Value Coupulas. A MSG is a copula C that is obtained as follows. 
Gonsider a multivariate distribution for a random vector [Xi, . . . , X„] having C = Cx 
as copula, take k iid copies of the vector given by 



i = 1, 



k 



and define 



Mkj = max 



Gonsider the copula 



for j < n. 



Cm{u) — CMk,i,Mk,2,-,Mk,S'^l^ . . .,Un) - C x^,X2,...,X„{u\, . . . , U^riY^^ 



The copula Cx is said to be Max-Stable if 



CxM = Cm(u 



namely if 



C{u'„...,u'^Y/' = C{u,,...,Un) 



which is precisely our definition of self chaining copula. 

One can see that the max stable copulas (or, equivalently, extreme value copulas, 
see Theorem 3.3.5 in [31]) follow the same characterization of our self-chaining copulas 
above, and it is known that the only archimedean max-stable copula is the Gumbel- 
Hougaard copula, see |TT] or Theorem 4.5.2 in [31J. 

Besides the other case of the Marshall-Olkin copula, featuring a discrete compo- 
nent, we may fully characterize Extreme Value (and hence Self-Ghaining) copulas by 
means of the Pickands function characterization. 
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Indeed, taking for simplicity again the bivariate case, basically our earlier ho- 
mogeneity characterization translates in the two-variables homogeneous function be- 
ing expressed as a one-variable function of the ratio. Such function consitutes the 
Pickands function of the copula, and Extreme Value Copulas (and hence Self-Chaining 
Copulas) can be fully characterized in terms of such function. For the multivariate 
case, the Pickands function is the restriction of the ra- variables homogeneous function 
to the unit simplex. For the details and the Pickands function properties see for 
example [13] or [26] . 

9 Conclusions, practical applications and further 
research 

This work investigates dependence across marginally exponentialy distributed inter- 
arrival times, for example default times in financial modeling or inter-failure times in 
reliability theory. 

In many applications one needs to simulate the random times progressively, check- 
ing "survival" in subsequent intervals 

[0,T), [T,2T), [2T,3T),..., [{N-1)T,NT] 

rather than at the final time NT directly. If one assumes exponential marginals for 
the random times and a given dependence structure (copula) in checking whether the 
arrival times "survive" each subinterval 

[{i — l)T, iT) conditional on having survived {i — l)T, i = 1, . . . , N, can one deduce the 
properties the dependence structure needs to feature to be consistent in the interval 
[0, NT] as a whole? As we have shown above, this condition is not easily satisfied. 
For example, iterating a Gaussian copula is not consistent with a Gaussian copula 
for the total [0, NT] interval, and in fact the inconsistent dependence obtained via 
iteration will be weaker than the dependence implied by the same Gaussian copula 
apphed directly and consistently to the single [0, NT], as we have seen in Example 
15.11 above. 

The self-chaining condition we introduce in this paper is necessary and sufficient 
to guarantee that a multivariate lack of memory property holds consistently with a 
given dependence structure. The dependence structure is consistent for the iterated 
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and terminal sampling. This means that we may choose the copula in the self-chaining 
family and then later decide whether we wish to sample times occurrences in a step 
by step fashion or in a one-shot fashion with the same copula without introducing 
inconsistencies. 

We point out that the self-chaining condition is the same as the extreme value 
copula condition and that, in the context of the Archimedean copulas associated with 
frailty distributions, the self-chaining condition forces the copula to be the Gumbel- 
Hougaard copula, thus providing a full characterization. The other known explicit 
case satisfying the self-chaining condition is outside the Archimedean family and is 
given by the Marshall- Olkin copula, featuring a discrete component and a positive 
probability that ti = T2. More generally, self-chaining copulas are characterized via 
a homogeneity condition for the log-copula, or a related partial differential equation, 
which translates into the Pickands function characterisation. 

Further research should focus on multivariate self-chaining copulas of dimension 
larger than 2 and with a flexible and rich parametric structure. It would be ideal if the 
parametric dependence were expressed via a matrix. Two such matrix-parametrized 
copulas are the Gaussian and Student-t copulas which, however are not self-chaining. 
A possible research idea is trying to extend the bivariate self-chaining limit of the t 
copula (see for example [13]) to the multivariate case. 
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